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Abstract 

We present a model which can explain semiquantitatively a number of the 
unusual properties of Yl>4 AS3 . The structural phase transition at T c ~ 300 K 
is described by a band Jahn- Teller effect of correlated electrons and is inter- 
preted as a charge ordering of the Yb ions. The low carrier concentration in 
the low-temperature phase follows from the strong electron correlations of the 
4f-holes on the Yb sites and can be viewed as self-doping of charge-ordered 
chains. The observed heavy-fermion behaviour is on a scale of T* ~ 50 K 
and is due to spinon-like excitations in the Yb 3+ -chains. The appearance 
of a second low-energy scale around 0.2 K is due to the Fermi energy of the 
low-density carriers. 
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The rare-earth pnictide Yb4As3 is a material with unusual physical properties 
At T c ~ 300 K, a structural phase transition has been observed which is accompanied by a 
charge ordering of the Yb ions and by discontinuities of the resistivity and the Hall coefficient. 
The high-temperature phase has the anti-Th 3 P 4 structure with space group I43d, while 
the one of the trigonal low-temperature phase is labelled R3c. Mofibauer spectroscopy 
data on 170 Yb give clear evidence that below 50 K the system contains about 20 % Yb 3+ 
ions characterized by a 4f 13 -configuration and a J = 7/2 ground-state multiplett. This 
fraction is close to the expected value of 25 % if perfect charge ordering takes place with 
the remaining 75 % Yb-ions in a Yb 2+ - or 4f 14 -configuration. No indication of magnetic 
ordering is found down to 0.045 K Instead, the material shows typical heavy-fermion 
behaviour with a linear specific-heat coefficient 7 ~ 200mJ/molK 2 and a correspondingly 
large spin susceptibility xs- The Sommerfeld- Wilson ratio is found to be of order unity, 
clearly indicating Fermi-liquid behaviour. However, below 2 K the susceptibility rises again, 
pointing to the presence of another low-energy scale . In agreement with the Fermi-liquid 
behaviour the resistivity is found to be of the form p(T) = p + AT 2 at low temperatures 
with a ratio A/^ u {y ~ 2) similar to that of other heavy-fermion systems 0. 

An unusual feature of Yb4As3 is that the heavy-fermion properties arise even though 
the density of charge carriers is extremely low. The low-temperature Hall coefficient i?n is 
positive, implying hole conduction and has a value of (eaRn) -1 = 7- 10 18 cm -3 . Thus the 
density of charge carrier is of order 0.001/Yb-atom. Although Yb4As3 is often referred to as 
a "low-carrier Kondo system" , it is very suggestive from the above that the heavy-fermion 
like behaviour of the system must be of different origin than in Kondo-lattice systems. (For 
reviews see, e.g., In fact, we are suggesting that it is unrelated to the Kondo effect. 

In the following we describe a model which can qualitatively (and semiquantitatively) 
explain a number of the unusual properties of the material, in particular the low carrier- 
density. It differs from the one suggested in which is based on a slightly depleted pnictogen 
band crossing a narrow Yb 4f band of width ~ 50 K. (The latter is assumed to have one hole 
per eight Yb ions.) Our model is based on a band Jahn- Teller effect of correlated electrons 
(cbjt) and rests on the following basic assumptions: 

• The fifty-six Yb 4f-bands have an overall width of ~ 0.2 eV which arises from effective 
f-f hopping via hybridization with pnictide valence states 0. In Yb4As3 all Yb-atoms 
are located on four families of chains oriented along the space diagonals of the cubic 
unit cell [|1(J. To reduce the complexity of the f-bands we make the assumption that 
they can be described by four degenerate one- dimensional bands associated with the 
chains. Such a model is rather similar to the Labbe-Friedel model for 3d states in A 15 
compounds where one has chains parallel to the cubic axes. Although such models 



may not be literally true due to interchain coupling, they describe the important aspect 
that a strain coupling to the one- dimensional band states may easily lead to distortions 
of the cubic structure and simultaneous repopulation among 4f band states of Yb4As 3 . 

There is a strong deformation potential coupling typical for mixed-valent systems 
which removes this degeneracy by a trigonal CBJT-distortion. When the Jahn- Teller 
transition takes place, the crystal shrinks in (lll)-direction, lifting the four-fold degen- 
eracy. The four equivalent chains are subdivided into one along the (lll)-direction, for 
simplicity referred to as short chain, and the three remaining chains referred to as long 
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chains. Since the Yb 3+ ions have a smaller radius than the Yb 2+ ones, it is natural 
to think of the cbjt transition being related to the ordering of Yb 3+ ions along the 
space diagonal (111). 

We describe the cbjt transition by an effective Hamiltonian of the form 

# = £ 4JW + e r E + ^L/W r , (la) 

where /i labels the summation over the /if = 4 different chains, and (ij) denotes a summation 
over nearest-neighbour sites along one chain. The fj create f-holes with spin a at site % of 
chain \i. We choose a bandwidth At = 0.2 eV as obtained from lda calculations H and, for 



simplicity, an effective spin degeneracy of 25 + 1 = 2. The second term in ( |Ta ) describes the 
volume-conserving coupling of the trigonal strain er (r = T 5 ) to the f-bands characterized 
by a deformation potential 

A, = A(^i-(l-<W^)- (lb) 

We assume that changes of the bandwidths due to the distortion are small and can be 
neglected. The third term in flTa]) is the elastic energy associated with the distortion. Here 
A^l is the number of sites in a chain and Co is the background elastic constant for one chain 
for which we choose cq/Q = 10 11 erg/cm 3 , where Q is the volume of a unit cell. (The lattice 
constant is a = 8.789 A.) Note that the bulk elastic constant c r = Acq. An elastic constant 
of this size is common for rare-earth systems, but for Yb4As3 it has not yet been measured. 
The Hamiltonian (1) does not yet contain the strong Coulomb interactions of holes at an 
Yb site. Therefore it is reasonable only above T c where the number of holes per Yb site is 
1/4 and somewhat below T c . At low temperatures the Coulomb interactions and the strong 
correlations which they imply are crucial (see below). However, for modeling the cbjt 
phase transition the Hamiltonian (1) is sufficient. The condition for a phase transition to 
occur is A 2 / (Atco) > 3. We choose A = 5 eV, which corresponds to a reasonable Gruneisen 
parameter Q = A /(At) = 25. From this, a transition temperature of T c ~ 250 K is obtained 
which is close to the observed T c exp ~ 300 K. 

Below the phase transition, one obtains a strain order parameter er(T) whose tempera- 
ture variation together with the band occupation n^(T) is shown in Fig. [I]. At T c , the four 
degenerate bands split with an associated shift of holes from the upper threefold degener- 
ate bands into the lower fourth band. The energy difference between the respective band 
centers is ||erA|. The equilibrium strain at low temperature is er = — A/(8c ) ~ —0.02. 
The increase in the hole occupation number of the lower band shows a similar behaviour as 
|er(T)| until at low T the upper bands are empty (of holes) while the lower band is becoming 
half-filled. 

As the lower band is approaching half-filling, the strong on-site f electron correlations 
become more and more important. Together with the coupling between short and long 
chains, this has important consequences as we will demonstrate. We incorporate the strong 
correlations by an on-site Hubbard interaction U which can be eliminated by going over to 
a t-J Hamiltonian of the form 
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FIG. 1. Temperature dependence of the strain order parameter er(T). Inset: Occupation 
numbers n^(T) of the four f-bands. 



# (1) = -t E rfLf^r + je (s! 1} ■ s« - i„w„f ) , (2) 



where J = At 2 /U and = E a /3 fla^apf^p, = E CT 4?- With [/ = 10 eV, we obtain 
J = 1 • 10 -3 eV. The projector V serves for the projection onto singly occupied sites. We 



make use of the so-called slave-boson approximation [|12] , |13| , representing the projector V 
by auxiliary bosons bj. The projection back onto the physical states of the now enlarged 
Hilbert space is done by means of Lagrange multipliers i\j for each site j. We introduce two 
Hubbard-Stratonovic fields Xij — (So- niafjio) an -d r j — (bj) to decouple the Hamiltonian 
(0), assuming that we work in a range of parameters where we can ignore Cooper pairing 
by setting (/m/jii — /iii/jit) — 0. ((S^ ) = as well.) Performing the uniform- saddle 
point approximation, Xij = X-> r j = r > = \ the Hamiltonian becomes after Fourier 
transformation 



E&AJW + N L g J X 2 - A (l - r 2 ) - \j (l - r 2 ) 2 ) , (3a) 

ka 



£ k = -2 [r 2 t + -^Jx) cos fc + A + ~ J (3b) 

The variational equations for the mean-field parameters lead to r 2 = 5 (where S is the 
deviation from half-filling of chain 1), A = 2t%, and x(T — 0) = ^ sin (| (1 — <5)). We observe 
that for low temperatures, the thermodynamics is essentially governed by a renormalization 
of the effective mass of the lower band, 

— = - * 3 T ^1.0-10 2 , (4) 

where is the bare band mass. This value for m* corresponds to a spinon bandwidth of 
25 K, which is close to the experimental T*. 
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This brings us to an important point, namely the determination of the quantity 5 in the 
limit of T — > 0. Without interchain coupling, the upper three bands are empty (of holes) 
in that limit, while the lower band is half-filled. The system should be a Mott-Hubbard 
insulator in that case. This changes when the interchain coupling is taken into account. In 
order to appreciate this point we neglect first the correlations. We introduce transversal 
hopping matrix elements t' ± er induced by the distortion between the short chain and the 
other three chains. These are again effective matrix elements, because the Yb-Yb hopping is 
viewed as taking place via intermediate pnictogen states. We have scaled t' ± with the strain 
to be able to restore the cubic symmetry at T > T c . Then we have to find new single-particle 
eigenstates <p\ Xa = E^lfif^AV an d the corresponding eigenvalues. The latter result is an 
additional repulsion of the lower and upper bands by an amount of 9t'^/(4A). The lower 
band is still half-filled, but when the strong correlations are taken into account, this does 
no longer imply a Mott-Hubbard insulating behaviour. Instead, the hole count on the short 
chain is reduced to 




^EI^WI 2 (4i CT M-i-^(^) , (5) 



since the single-particle eigenstates <pk\a contain contributions oc t' ± /A from the long chains. 
The Hubbard-f/ has an effect only in the short chains where the hole concentration is nearly 
1, but not in the long chains where the hole density is very small. Therefore the system 

acts like being self-doped with a doping parameter 5 = || . Taking 5 ~ 0.6 ■ 10~ 3 from 
experiment, we have t' ± ~ 0.1 eV. Note that the self-doping found here differs from the one, 
e.g., in YBa2Cu3C>7 where holes are transferred from the Cu 3+ ions in the chains to the Cu 2+ 
ions in the planes because several bands are crossing the Fermi level. 

With the picture outlined above one can understand, at least qualitatively, a number 
of properties of Yb4As3. The large 7-coefficient results from the (spinon-like) excitations 
within the short chains. We obtain it formally from the large mass enhancement @. The 
same holds true for the spin susceptibility. At this point it should be noted that our model 
contains a second temperature scale, which is given by the Fermi energy Ep of the holes 
in the long chains. For 5 ~ 0.6 ■ 10 _3 /Yb-atom we find ~ 0.2 K. This might explain 
the low-temperature behaviour of xs{T), because the holes in the long chains give rise to a 
Curie-like susceptibility except at k^T <C when the Fermi gas becomes degenerate. 

As far as transport is concerned, we expect that the same heavy quasiparticles responsible 
for the large 7-coefficient are also governing the low-T behaviour of the resistivity. The 
resistivity of the short chains is presumably higher than the one of the holes in the long 
chains. In a crystal with different domains the one of the short chains will therefore make 
the largest contribution to the measured resistivity. This might explain why the prefactor 
A of the T 2 -term leads to a value for the ratio A/7 2 like in other heavy-fermion systems 



14| . As far as p(T) for T < T c is concerned, it must increase with decreasing temperature 
because of the depletion of holes in the upper bands. Thus p(T) must have a maximum in 
order to match on to the po + AT 2 behaviour at low T. Concerning the Hall coefficient, 
we expect that in the presence of different channels with different conductivities the one 



with the largest a makes the biggest contribution [15] . This would imply that the holes 



in the long chains are most effective here. Since the hole density increases with rising T, 
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the Hall coefficient decreases accordingly. The initial increase for small temperatures is 
possibly an effect beyond the simple one-band model expression for the Hall coefficient. It 
is worth noticing that in the high-temperature phase the experimental Hall coefficient is 
(eaRn) -1 — 1 • 10 21 cm~ 3 . With a volume Q = 6.8 • 10~ 22 cm 3 per unit cell, this would 
correspond to 0.25 holes per formula unit Yb4As3 provided we have a one-band system. 
With four independent bands the application of a generalized theory for the Hall coefficient 
n| yields one hole per formula unit, as expected from a chemical approach using valence 



electron counting. 

In conclusion we may state that the model presented here for Yb4As3 explains semiquan- 
titatively the most important experimental findings in a natural way. Among them are the 
low carrier concentration in the low-temperature phase, the heavy-fermion behaviour with 
an effective mass ratio m*/m ~ 100 as well as the appearance of a second low-energy scale 
of order 0.2 K. 
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